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The IceCube observatory is the first cubic kilometre scale instrument in the field of high-energy neutrino astronomy and cosmic 
rays. In 2009, following five successful deployment seasons, IceCube consisted of 59 strings of optical modules in the South Pole 
ice, together with 118 air shower detectors in the IceTop surface array. The range of physics topics includes neutrino signals from 
astrophysical sources, dark matter, exotic particle physics, cosmic rays, and atmospheric neutrinos. The current IceCube status and 
selected results are described. Anticipated future developments are also discussed, in particular the Deep Core low energy subarray 
which was recently deployed. 
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1. Introduction 

At the IceCube observatory, the clear South Pole glacial ice 
is exploited to detect Cherenkov light emitted by charged par- 
ticles created in neutrino interactions. An array of digital op- 
tical modules (DOMs) equipped with photomultipliers extends 
throughout a volume of approximately 1 km 3 , with the central 
region more densely intrumented. Scheduled for completion 
in 2011, IceCube allows the detection of neutrino interactions 
above a few tens of GeV, with a maximum sensitivity in the 
TeV-PeV range. 

Extraterrestrial neutrinos are expected from a variety of 
sources. Generically, high energy cosmic ray particles will 
give rise to neutrinos in their interactions with matter or radia- 
tion, either in the immediate vicinity of the acceleration sites or 
while traveling through space. It is thus possible to estimate the 
neutrino flux that corresponds to the observed flux of cosmic 
rays Q. Such considerations have dictated sensitivity require- 
ments which translate into a cubic kilometre detector. 

Owing to their low interaction cross-section and lack of elec- 
tric charge, neutrinos propagate over vast distances, quite undis- 
turbed by matter, radiation, and magnetic fields. This is in 
contrast to photons, which can be absorbed by matter, and 
which cannot travel cosmological distances if their energy is 
above the TeV range, because of absorption by background ra- 
diation. Similarly, protons of the highest energies can be ab- 
sorbed through the GZK mechanism [2], while at lower ener- 
gies they are deflected by magnetic fields and do not point back 
to the source. It is therefore possible that there are "hidden" 
sources which can only be detected via their neutrino emission. 
In cases where emission of gamma rays or high energy cos- 
mic rays is seen, a detection of neutrinos could elucidate the 
nature and location of cosmic ray accelerators. (Gamma rays 
from 7T° decay would be accompanied by neutrinos from n + jn~ 
decay.) Promising candidate acceleration sites are relativistic 
ejecta from active galactic nuclei (AGNs) or gamma-ray bursts 



(GRBs). There are also objects within our galaxy, such as su- 
pernova remnants, pulsar winds, or microquasars, in which ac- 
celeration processes could yield neutrino emission. 

Even if individual neutrino sources are not seen, their com- 
bined diffuse flux may stand out above the atmospheric neutrino 
flux. The study of the atmospheric neutrino flux is therefore an 
important topic. It provides a useful calibration point and may 
in addition yield results relevant to the intrinsic properties of 
neutrinos. 

There are abundant indications of dark matter in the 
form of non-relativistic weakly interacting massive particles 
(WIMPs) 0. These could constitute a class of "hidden" 
sources, detectable via neutrinos, as they annihilate in the centre 
of the Sun or other dense regions (4j|3. Promising candidate 
WIMPs are the neutralinos of supersymmetric models or 
Kaluza- Klein excitations in models of extra dimensions (6). 

Section [2] describes the IceCube detector and its perfor- 
mance, and selected analysis results are discussed in section 
O Future enhancements are discussed in section [4] 

2. The IceCube detector 

Figure [T] illustrates the layout of the IceCube observatory in 
the polar icesheet. Also shown is the AMANDA array Q, 
which continued to operate as a part of IceCube until it was 
finally closed down in May 2009. In addition to the deep ice 
detector, there is a surface array, IceTop, for detecting cosmic 
ray air showers. Each IceTop station is located above a string, 
and consists of two tanks. Each tank is filled with clear ice and 
monitored by two DOMs running at different gains to achieve 
the required dynamic range. 

IceCube construction takes place during the austral summer. 
A 5 MW hot- water drilling system is used to melt holes in the 
ice down to a depth of 2,45 km. In each hole a string consisting 
of a cable carrying 60 optical modules is deployed, and the wa- 
ter in the hole then re-freezes. During the recent deployment 
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Figure 1 : The IceCube detector. Filled circles indicate the 59 strings completed 
for the 2009 data taking season. The shaded region near the centre is the Deep 
Core subarray, and the one to the right shows the position of AMANDA. 

season, the time between completed deployments was about 
two days. When completed in 2011, IceCube will consist of at 
least 80 strings and 80 IceTop stations. Fifty-nine strings were 
in operation during 2009, and an additional 20 were deployed 
in the summer season of 2009/2010. 

Figure [2] shows an IceCube optical module (8). The light- 
sensitive element is a 25 cm photomultiplier tube (PMT), look- 
ing down. This is optically coupled to a surrounding pressure 
sphere via transparent gel. The sphere also contains a high- 
voltage generator, 12 LED flashers for calibration purposes, 
communication and time-calibration electronics, and two dif- 
ferent types of digitization hardware. Three Analogue Transient 
Waveform Digitizer (ATWD) channels with different gains are 
used to digitize the PMT output with nanosecond precision dur- 
ing 400 ns. There is also a fast ADC running at 40 MHz during 
6,4 jis. The DOM is triggered by a discriminator coupled di- 
rectly to the PMT output, and the input to the digitizers is de- 
layed by 75 ns to allow the capture of the complete waveform. 

An important component in the IceCube detector is the ice 
itself. Its optical properties have been measured in situ using 
light emitting devices |9], and the ice model and its implemen- 
tation in light propagation and reconstruction software are be- 
ing gradually refined. Scattering dominates over absorption, 
with an overall attenuation length in excess of 20 m in the up- 
per regions of the detector, and about twice that value below a 
depth of 2100 m. In between these two regions there is a dusty 
layer with strong scattering and absorption. 

IceCube is designed to detect neutrinos of all flavours and 
distinguish extraterrestrial neutrino signals from the back- 
ground. The remainder of this section contains a discussion 




Glass Pressure Housing 

Figure 2: A schematic drawing of an IceCube optical module (DOM). Figure 
from reference [8] 

of the principles involved, as well as some examples of detector 
performance. In all cases, detection relies on the emission of 
Cherenkov light, at an angle of 41° from the direction of prop- 
agation of charged relativistic particles. 

When a muon from a charged current (CC) interaction 
passes through the instrumented volume it gives rise to a char- 
acteristic "track" consisting of large pulses in DOMs close to 
the muon trajectory, separated in time by the muon flight time. 
Its direction can therefore be well determined. The point spread 
function for well reconstructed muon neutrino events is shown 
in Figure [3] Sub-degree resolution has been achieved already 
with the 40-string configuration. Muons which pass close to 
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Figure 3: The cumulative point spread function for well reconstructed events, 
including the angle between the incoming neutrino and the muon. The improve- 
ment between the 22- and 40 string case is not only due to the larger detector, 
but also to an improved reconstruction algorithm. Figure from reference fTol 

the edge of the instrumented volume, or outside, can also be 
detected, especially at high energies. In this case, however, the 
event reconstruction and background rejection is more difficult, 
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and relies more heavily on the description of the emission of 
Cherenkov light and its propagation through the ice. 

To achieve sensitivity to extraterrestrial neutrino sources, one 
must also eliminate the overwhelming background of downgo- 
ing muons from cosmic ray air showers. The traditional ap- 
proach is to use the Earth as a filter and select up-going events. 
There will still be a large background from downgoing muons 
which are mis-reconstructed as upgoing. By requiring well re- 
constructed events, however, one can define a clean sample of 
up-going neutrinos, as illustrated in Figure [4] 
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Figure 4: Muon trigger rates in the 22-string IceCube configuration, as a func- 
tion of the cosine of the zenith angle. The rate in data (red squares) is compared 
to a simulation including both atmospheric muons and muons from interacting 
atmospheric neutrinos. Requirements to select well reconstructed events were 
applied, and the rates are corrected for the efficiency of these. Figure from 
reference fTTI 
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Figure 5: Reconstructed vs true logarithm of the muon energy, at the point of 
closest approach to the centre-of-gravity of hit modules. Note the logarithmic 
intensity scale. 
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The energy spectrum of atmospheric neutrinos falls rapidly, 
as ~ £"" 3 ' 7 , whereas neutrinos from atrophy sical sources are 
expected to have a much harder spectrum (~ E~ 2 ). An en- 
ergy estimator will therefore provide additional discrimination 
in the search for extraterrestrial point sources. Furthermore, it 
is essential in the search for a diffuse signal of extraterrestrial 
muon neutrinos. For muons above a few TeV, the amount of 
Cherenkov light scales roughly linearly with energy, and the 
total number of triggered DOMs has been used as an energy 
proxy. Recently, a more direct method for measuring the energy 
loss along the muon trajectory has been developed (T71 [T8lL 
with a resolution in log 10 (£^/GeV) of around 30% in the range 
between 1 TeV and 100 PeV, see Figure [5] 

Two recent observations of downgoing atmospheric muons 
illustrate that the reconstruction methods work. One is the de- 
tection of the reduction of the cosmic ray flux caused by ab- 
sorption in the Moon lfT2l . see Figure [6] Systematic biases in 
the reconstruction could affect the position of the minimum. 
No such effect is seen in the analysed data. Increased statis- 
tics will render the test more stringent. The other result is the 
observation of a global anisotropy in the cosmic ray arrival di- 
rections tl3\ . This is the first such measurement in the Southern 
Hemisphere, and it agrees well with measurements done in the 
North d El 

Neutral current interactions give cascade events in which no 
particle travels a distance comparable to the spacing between 



Figure 6: Angular distance from the Moon for events in a 1,25° declination 
band centered on the Moon. A 5,2cr deficit is seen in the central bin. The figure 
is based on eight months of data taken with the 40-string configuration in 2008. 



Bartol + Naumov RQPM 
AMANDA-II 2000-6 v,, 13B7 d 
Razzaque GRB Progenitor 2008 
Blazars Stecker 2005 
Baikal 1038 d cascades 
AMANDA II 3 year v u x 3 
IC22 Vj,x3 ICRC2009 Sensitivity 
Waxman Bahcall 1998 x 3/2 

' +v e , 2001 
HiResv t +v e x3/2, (2008) 
IC22 EHE(all flavor),242 d ICRC2009 



Honda + Sarcevic Min 
IC22 Atmo. ICRC2009 
Waxman Bahcall Prompt GRB 
BL LACs Mucke et all 2003 
AMANDA-II 2000-4 cascades 
AMANDA/TWR 3 year UHE (all flavors) 
IC40 v„x3 Preliminary Sens. 
Rice (all flavors),(2006) 854 d 
PAO v,x3,2 yr ICRC 2009 
ANITA-2008 (all flavors), 35 d 




log io E, [GeV] 



Figure 7: Diffuse neutrino flux predictions, measurements, and limits. All 
results are scaled to apply to the total neutrino flux, assuming that the fluxes are 
equal for all three flavours. See text for details. 
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detector modules. This is also true for v e CC events below the 
EeV range, and for v T CC events below a few PeV. For such cas- 
cades, the light is emitted from a small region, but preferentially 
at the Cherenkov angle of 41° with respect to the shower axis. 
Making use of this information, and detailed knowledge about 
light propagation in the ice, an angular resolution of 30 to 35 
degrees is expected for cascades in the energy range between 
10 TeV and 10 PeV ESI. The resolution in log 10 (£/GeV) is 
~ 13% for cascades contained in the instrumented volume fT6ll . 
As the array grows, the volume where a cascade can be identi- 
fied increases faster than the detector volume, and much faster 
than the area exposed to a throughgoing muon. Therefore the 
cascade channel, with its better energy resolution, is becoming 
more important as the detector grows. 

At higher energies electromagnetic showers initiated by elec- 
trons from y e interactions become elongated owing to the LPM 
effect (191 [20), which makes a better angular reconstruction 
possible I2H . For v T energies in the PeV range and above, the 
distance travelled by the r from a CC interaction leads to sev- 
eral signatures [22]. These are of special interest as there is no 
terrestrial v T background flux, whereas astrophysical neutrino 
fluxes are expected to contain a mix of all flavours. 

3. Selected results 

Figure [7] shows two preliminary results on the diffuse neu- 
trino flux, obtained from the IceCube data taken in 2007 with 22 
strings (IC22). The triangles with error bars show a preliminary 
energy spectrum, obtained using an unfolding method ifTTl . 
Agreement is found with previous AMANDA unfolding results, 
and with the calculated atmospheric neutrino flux. The con- 
nected filled circles show a quasi-differential limiQat 90% con- 
fidence level, obtained in a preliminary low background search 
for Extremely High Energy (EHE) events [23]. No candidates 
were found. 

Also shown are projected sensitivities of the searches with 
IC22 and IC40 (the 2008 40-string configuration) for extrater- 
restrial diffuse fluxes falling as E~ 2 fISl . They are drawn as 
constant lines with a length corresponding to a region which 
would contain 90% of the selected signal events. The figure 
also shows results obtained with AMANDA. These are an un- 
folded energy spectrum for atmospheric neutrinos, based on 
seven years of muon data (24), limits on a diffuse E~ 2 flux ob- 
tained using muons [ 25 ] and cascades [26], and a limit on the 
flux of Ultra High Energy (UHE) neutrinos based on three years 
of AMANDA data (27). 

Figure [7] also includes some results from other detectors. 
They are, in order of increasing energy, the Lake Baikal neu- 
trino observatory, the RICE South Pole radio Cherenkov ar- 
ray, the Pierre Auger Observatory (PAO), the HiRes air fluo- 
rescence experiment, and the ANITA radio Cherenkov detec- 
tion antenna. Also shown are the Waxman-Bahcall flux derived 



lr This is the limit on a spectrum which falls as E 2 over one decade in en- 
ergy, and is zero outside this range. 
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Figure 8: Significance map showing - log 10 p for points in the northern sky, 
based on a search for point sources in seven years of AMANDA data. 

from the observed cosmic ray spectrum [ 1 ] and calculated con- 
tributions to the diffuse neutrino fluxes from atmospheric neu- 
trinos, gamma-ray bursts (28), active galactic nuclei |29lL and 
the GZK process l30l . 

Some of the neutrino sources contributing to a diffuse flux 
may be individually seen as point sources. Since the energy 
spectrum of the signal is expected to be harder than that of the 
atmospheric neutrino background, a likelihood ratio hypothesis 
test including energy and direction is the natural approach |3T1 . 
Repeating this for a very fine grid of source directions gives a 
map of the /7-value (the a priori probability, in the absence of a 
source, to obtain a result which is at least as signal-like as the 
one actually obtained). 

Figure [8] shows a sky map obtained from seven years of 
AMANDA observations [32]. The maximum significance is at 
S = 54°, a = 171°, where the /7-value is p = 7,4 • 10" 4 . The 
probability of obtaining such a /7-value or lower somewhere 
in the sky, determined by randomising the right ascension, is 
p = 95%. Hence there is no indication of a point source signal. 

The point source search in the data taken with the IceCube 
22-string configuration lf33l was about twice as sensitive as 
the AMANDA one, and showed a conspicuous deviation from 
background with p = 7 • 10" 7 at 6 = 11,4°, a = 153,4°. How- 
ever, the final /?-value, accounting for the multiple trials in scan- 
ning the sky, is p = 1,3%, which is a very marginal significance. 

A point source analysis of six months of data taken with 
40 strings (out of 13 months in total) ifTUl yielded the sky 
map shown in Figure [9] Here the most significant point is at 
S = 15,4°, a = 115,0° and the p- value, accounting for mul- 
tiple trials, is 61%. The map includes the southern sky, where 
the background of downgoing muons completely dominates the 
atmospheric neutrinos. A zenith dependent cut on the muon en- 
ergy parameter was imposed to ensure a nearly uniform zenith 
distribution of the background. As a result, the search in the 
southern hemisphere is sensitive primarily to neutrino energies 
in the PeV range, as compared to TeV energies for the northern 
hemisphere. The better angular resolution for higher energy 
can be clearly seen in the figure. A search above the horizon 
has also been performed using the 22-string data (34). 

In each of the point source analyses mentioned above, source 
candidates on a pre-defined list have been tested, in addition to 
the very fine grid search. The most significant result is for the 
Geminga pulsar in the seven-year AMANDA data (p - 0,9%). 
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Figure 9: Significance map for a point source search in six months of data 
taken with the IceCube 40-string configuration. There are 6796 events from 
the northern celestial hemisphere, and 10981 from the southern one, where the 
atmospheric muon background dominates. 
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Figure 10: Calculated prompt neutrino fluences for 41 GRBs searched for 
neutrino emission in the 22-string data. Also shown is the sum over all bursts 
(thick solid line), the Waxman-Bahcall spectrum | 42 1 for a single burst (thick 
dotted) and for 41 bursts (thick dashed). 



For the IceCube searches, however, Geminga gives downward 
fluctuations (p > 0,5). The point showing the excess in the 22- 
string data was included as a candidate in the 40-string search, 
where it showed a downward fluctuation. 

A search for time clustering of events from candidate sources 
has also been carried out using the 22-string data [ 35 ], and there 
have been efforts to correlate IceCube observations in time with 
flares observed in the electromagnetic spectrum (36). There is 
also a programme to provide online triggers for air Cherenkov 
telescopes |37], and for fast optical follow-up observations of 
transient neutrino emission, as expected from e.g. supernovae 
or gamma-ray bursts (38). 

In addition, the recorded data are searched for neutrino 
events coincident in time and direction with GRBs detected by 
satellites (39j|40l|4l]]. The approach has recently been refined to 
model the emission of individual bursts based on the observed 

Limits obtained (39) are 



parameters, as shown in Figure 10 
two orders of magnitude above the fluence sum in Figure [TO 
and about a factor five above a predicted lower energy precur- 
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Figure 1 1 : Limits on the spin-dependent neutralino cross-section obtained with 
IceCube and AMANDA compared to results from Super-Kamiokande and di- 
rect detection experiments (cited in reference 1431 . Also shown is the region 
corresponding to allowed supersymmetric models found using DarkSusy H51 
for current limits on the spin-independent cross-section, and the region which 
would still be allowed if these limits improved by a factor 1000. 



sor fluence|28]. A study was also made l39l of the expected 
sensitivity of the full IceCube detector, with future GRB sam- 
ples from the Fermi and Swift satellites. It indicates that flu- 
ences, predicted from the assumption that GRBs are the domi- 
nant sources of ultra-high energy cosmic rays, will be detected 
or ruled out after some years of running IceCube l39l . 

In conclusion, no evidence for steady or transient neutrino 
point sources has been found so far. 

The indirect detection of dark matter via neutrinos from 
WIMP annihilations is another important aim of IceCube. The 
nearest candidate sites for gravitational trapping are the centres 
of the Earth and of the Sun. Since the light elements dominate 
the Sun, the spin-dependent scattering cross-section is more rel- 
evant for WIMP capture in the Sun than in the Earth. For the 
same reason, searches for an indirect signal from capture in the 
Sun complement those in direct detection experiments, which 
are primarily sensitive to the spin-independent cross-section. 

The 22-string IceCube data taken with the Sun below the 
horizon have been analysed, using a likelihood ratio method in- 
volving the angle relative to the Sun, to test hypotheses of neu- 
tralino annihilations l43l . Neutralino masses in the range from 
250GeV/c 2 to 5TeV/c 2 were selected, and two extreme anni- 
hilation channels were investigated, the "hard" (W + W") and the 
"soft" (bb). The data show a slight deficit in the direction of the 
Sun |43l . Following the procedure of referencel44l the resulting 
limits on the muon flux in the ice were converted into the limits 

The 



on the spin-dependent cross-section shown in Figure 1 1 
figure also shows limits obtained using AMANDA data |46l . 
One set was obtained using the data sample and techniques 
from the seven-year point source search l32l . and another from 
a dedicated analysis of three years of data 1471 . It can also 
be seen from Figure [TT] that the limits on the spin-independent 
cross-section from direct detection searches do not appreciably 
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Figure 12: Limits on the spin-dependent cross-section of the lightest Kaluza- 
Klein particle (LKP) and the region allowed in minimal models of extra dimen- 
sions. A (i) is the relative mass splitting between the first quark excitation and 
the LKP. The dark narrow (blue) region is obtained when consistency with the 
relic density inferred from WMAP measurements 1491 is required within one 
sigma. 



constrain the supersymmetric model space in this projection. 

The negative result of the search in the IceCube 22- string 
data was also interpreted [ 48 ] in the framework of Kaluza- Klein 
theories with universal extra dimensions, with the WIMP anni- 
hilating in the Sun being B (1) , the first excitation of the weak 



hypercharge gauge boson. Figure [12] shows the cross-section 
limit obtained, and the region allowed in minimial models. The 
limit from the seven-year AMANDA point source sample is al- 
most identical. 



4. Future 

There are several initiatives to benefit from the presence of 
the IceCube array and extend the range of physics topics which 
can be addressed. 

The IceCube volume is too small to collect a sizeable sam- 
ple of neutrino events at EHE energies, because of the low flux. 
Acoustic and radio emissions from neutrino interactions are at- 
tenuated much less than light, and could be used to increase the 
effective volume at high energies in a hybrid approach [50]. 

The South Pole Acoustic Test Setup [ 51] has been used to 
study sound propagation and the acoustic environment in the 
ice using pingers and receivers deployed in IceCube holes. The 
preliminary results indicate an attenuation length of around 
300 m, which is more than an order of magnitude below the 
predicted value (52). 

Within IceCube there has also been a successful development 
programme, AURA, to test methods for radio detection in the 
ice (53). A proposal for an ~ 80km 2 array (Askaryan Radio 
Array, ARA) has been submitted. In addition, the possibility of 
a surface radio array to detect air showers was studied (54). 

At the other end of the energy scale, WIMP annihilations 
would probably give neutrinos in the energy range from GeV to 
TeV. The original IceCube design is not optimal for these ener- 
gies. Additional funding was secured, however, for six strings 
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Figure 13: The effective area for muon neutrino detection at trigger level of the 
originally planned 80-string array, and with the addition of the Deep Core 



added at the centre of IceCube. These are equipped with PMTs 
with enhanced quantum efficiency, about 30% higher than that 
of standard DOMs. The modules are placed closer to each other 
than on standard strings (7m as compared to 17m). The instru- 
mented region is between 2100 m and 2450 m (below the main 
dust layer), where the exceptionally clear ice, together with the 
denser instrumentation, contributes to an overall ten-fold in- 
crease in photon collection efficiency. The lower portions of 
the central thirteen strings, including the aforementioned six, 
constitute the Deep Core subarray (55) , with a fiducial mass of 
about 13 MT. It is surrounded by at least three layers of instru- 
mentation on all sides. Simulations indicate that an algorithm 
rejecting events with early hits outside the Deep Core will be 
able to reduce the atmospheric muon background by a factor of 
10 6 , bringing it down to the level of atmospheric neutrinos (56) . 



Figure 13 shows the IceCube effective area with and without 
the additional Deep Core strings. The improvement is signifi- 



cant at low energies. The thick dotted line in Figure 11 shows 
the expected sensitivity in the solar WIMP search ten years af- 
ter the completion of IceCube, using data taken with the Sun 
below the horizon. The sensitivity below 100 GeV/c 2 is almost 
exclusively due to the Deep Core. 

The physics potential of the Deep Core goes well beyond 
the search for Dark Matter (57). The possibility to reject at- 
mospheric muons entering through the surrounding detector ar- 
ray will make it possible to study neutrino emission from the 
southern hemisphere without requiring PeV energies. Further- 
more, the atmospheric neutrino spectrum in the range between 
10 GeV and 1 TeV will be accessible, which will make it pos- 
sible to conduct high-statistics studies of neutrino oscillations 
at ~ 10 GeV over baselines comparable to the size of the Earth. 
It will be possible to study disappearance at higher energies 
than previously done, and possibly v T appearance. It has also 
been suggested that the Deep Core might have some sensitivity 
to the neutrino mass hierarchy (58). 

The additional six Deep Core strings were successfully in- 
stalled during the 2009/2010 deployment season together with 
14 standard IceCube strings, bringing the total to 79 strings. 
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